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Methods are presented for obtaining upper and lower pointwise hounds for the deflection (and its 
derivatives) in certain physically interesting elastic plate problems. 



1 . Introduction 

In a recent paper [1] ' the authors presented a new method for obtaining pointwise bounds 
in the first boundary value problem for elastic plates. The method is based on a mean value 
inequality and an a priori inequality for biharmonic functions. In this paper we use the same 
techniques to obtain pointwise bounds in terms of known data for the deflection and its derivatives 
in the second boundary value problem for an elastic plate, i.e., the displacement and bending 
moment prescribed on the boundary of the plate. In addition we derive pointwise bounds for the 
solution (and its derivatives) of the third boundary value problem, i.e., the bending moment and 
reaction prescribed on the plate boundary. In each of these problems bounds at an interior point 
for the deflection and any derivative up to and including the third are obtained in a straightforward 
manner. The inequalities obtained, bound an arbitrary (sufficiently smooth) function at a point 
in a region in terms of the data of the given problem. Thus if one applies the inequality to the 
difference between a desired solution (with known data) and an arbitrary trial function, it yields 
a bound for the error in the approximation of the unknown solution by the known trial function. 
The trial functions need only possess piecewise continuous fourth derivatives and are not required 
to satisfy any differential equation or boundary conditions. If the applied loading is sufficiently 
smooth in a neighborhood of the point at which bounds are sought, it is possible to obtain in addi- 
tion to the above mentioned results, bounds for higher derivatives of the deflection in either of the 
two problems considered here. We note the somewhat unexpected fact that it is considerably 
easier to obtain bounds in the third boundary value problem than in either of the first two. 

. In a subsequent paper the authors wih present methods for obtaining pointwise bounds in 
mixed boundary value problems for elastic plates. 

2. Notation and Definitions 

Throughout this paper we assume the plate to occupy a two dimensional bounded region R 
with boundary C. If R is not polygonal it is necessary (for our method) in the second boundary 
value problem to assume that C have bounded curvature. We make use of the summation con- 
vention throughout this paper and employ a comma to denote differentiation, e.g., W, i = dW/dx i . 
The symbol cr is used to denote Poisson's ratio and D to signify the plate rigidity. The symbol A 
denotes the Laplace operator and A 2 the biharmonic operator. Additional notation will be defined 
as need arises in the text. 



*University of Maryland, College Park, Md. Part time at National Bureau of Standards, Washington, D.C 
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145 



3. Second Boundary Value Problem for Elastic Plates 

The boundary value problem considered in this section may be stated as: 

in ti and 

u = g,M(u) = h (3.2) 



on C. The quantities/, g, h are prescribed data and M(u) is the bending moment; i.e., 

M(u) = -D 



Au-(l-cr)(^+ 



d 2 a ■ 1 du 
ds 2 p dn 



(3.3) 



In (3.3) ^-^ denotes the second tangential derivative of u, du/dn the normal derivative and p is the 
ds z 

radius of curvature on C. 

It is possible to obtain bounds for the strain energy, but since this quantity is not usually of 
physical interest we do not present a bound for it here. We indicate now a method for obtaining 
pointwise bounds for u and its derivatives. The result may be stated as follows: 

THEOREM I: For any biharmonic function V with piecewise continuous fourth derivatives in 
a domain R whose boundary C has bounded curvature, the following a priori bound holds for any 
pth order derivative of Y at a point O in R: 

|FW(0)| 2 ^Ki(p)£ V*ds + K 2 (p) i (j^ 2 ds + K 3 (p) j> (0) 2 ds + K 4 (p) & [M(V)] 2 ds, (3.4) 

where the constants /£i(p), . . ., K 4 (p) are explicitly determined. 

Before proving this theorem we indicate how such an a priori inequality leads to pointwise 
bounds for u and its derivatives in R. If we denote by <p any particular solution of (3.1) and set 
V—u — ip — B, where B is biharmonic, then by choosing the data of B to approximate that of u — <p 
(the usual Rayleigh-Ritz technique 2 ) we obtain pointwise bounds for u and its derivatives. 

In proving theorem I we make use of mean value inequalities derived in [1], i.e., if S r denotes 
the interior of a sphere contained in R and of radius r about O, then 

l^0)| 2 ^j^ V*dxdy*±fv*dxdy (3.5) 

\yii)(0)\*^2£- f VHxdy^^ [ V 2 dxdy (3.6) 

-rrr Js r trr 4 J R 

with similar expressions for the higher derivatives (see [1]). Thus in order to establish (3.4) it is 

sufficient to obtain a bound I V 2 dxdy in terms of the quantities on the right of (3.4). 

J R 

To this end we decompose V as 

V=B + H (3.7) 

where 

A//-0 (3.8) 

in R and 

H=V (3.9) 

on C. Then 

&B = (3.10) 

in R and 

B = 



(3.11) 



M(B)=M(V)-M(H) 



1 Cf., Handhuch der Physik, 1, Mathematische Methoden, Springer Berlin, p. 278 (1956). 
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on C. Clearly 



(3.12) 



[ l-dxdy^ 2i f B 2 dxdy+ f IPdxdyV 
We first derive a bound for I B 2 dxdy. To do this we introduce an auxiliary function ^satisfying 



ft and 



AW=B 
W=0,M(W)=0 
on C. Then using the divergence theorem we ohtain 

[ B 2 dxdy = ( BAWdxdy = -jjj> M(B) ^ds. 

Thus by Schwarz's inequality 



!,***«[»$< 



UmBm-idsltffi 



(3.13) 

(3.14) 

(3.15) 



(3.16) 



where / is any positive function on C. 

In order to hound the last integral in (3. 16) we introduce a vector field with components/' such 
that Jhn > on C. Then setting 

t=fm (3.17) 

we have 

2 I /'7F, /Arrf.rd>'=*/f^-Vr/.s+ (/,i:8{-2/,f)r,;r,^dy. (3.18) 

This is one form of an identity due to Rellich (ef. [5]). 

If the derivatives of/' are hounded in R, we ohtain by applying the arithmetic-geometric mean 
inequality to (3.18) 



if)*'- 



(3.19) 



■'A fw,iW 9i dxdy+a j(AW) 2 dxdy. 
Here c is just some bound for the largest eigenvalue of the matrix 2f\ ,•— ff k d\ i.e., 

f (2/, t-fStfQW. iW, idxdy *zc( V, ,JP, ,-dzdy, (3.20) 

and 

& = Max (/!/*). (3.21) 

In (3.19), a is any positive constant. For a star shaped boundary we may choose f i = x i in which 
case we may take c = 0. Now, since Jf 7 vanishes on C, it is clear that 



f W 2 dxdy ^7- f r, , JT, ;</.r</y s= -^ f 
J ft A.1 J/.' Aj J// 



dxdv^y, J (MVfdxdy 



(3.22) 



where Xi is the first eigenvalue of the fixed membrane problem for ft. 3 
The strain energy E{W \ W) is given by 



E(W 9 ^> = §{[ [<r(AW?dxdy+(\ -o)W,vW 9 v]dxdy\- 



(3.23) 



:! For a discussion of the fixed membrane problem see for instance Courant, R., and Hilbert, D., Methods <»f Mathematical Physics 1, Interscience (1953). 
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It follows easily, since W, U W, u ^ | (AW ) 2 , that 



jjE{W 9 W) 



From (3.22) 



[(AW) 2 dxdy^ T ^— 

JR 1+CT 

and, since W and M{W) vanish on C, -| £(JT, JP) = | WAWdxdy. Hence 

f (AWfdxdy ^ -r^— f WAWdxdy, 
Jr \L~rO') Jr 

I W 2 dxdy *£ 75 f (AWfdxdy. 

JR A.1 J« 

We thus obtain, using Schwarz's inequality on (3.25) and inserting (3.26), 

f (AW) 2 dxdy^ \ n J rr I' f (A 2 W) 2 dxdy. 

JR |_(H-0")AiJ Jfl 

Combining (3.19), (3.22) and (3.26) we now have (choosing a = [i/A.i] % ), 



i-m-a. 



-cr) 2 X3 



(c + 2 



[6X1] '^ f 



{t^WYdxdy. 



(3.24) 

(3.25) 
(3.26) 

(3.27) 
(3.28) 



In (3.28) we may use any lower bound for Xi, e.g., that obtained from the isoperimetric inequality 
of Faber [3] and Krahn [4] 

X, *£ njilA. (3.29) 



Here j is the first zero of the Bessel function J and A denotes the area of R. We now set 

4 



A, = - 



5=,(c + 2[6Xxl^) 



(l + o-) 2 X3 

where Xi is any lower bound for Xi. Then from (3.16) we obtain (recalling (3.13)) 

J B 2 dxdy^j±b [M(B)]*ds 

From (3.3), (3.9), and (3.11) we have 

M(B)=M(V)-D(l-a) 
Thus 



(3.30) 



(3.31) 



dW 1 dH 
ds 2 p dn 



L B2dxd ^ 2 M 



M(V)-D(l-a) , 



9V 
ds 



ds + D 2 {\-<T) 2 <\>±- 2 [™) ds 



Wd/ft 2 

c p 2 \dn , 



We assume that the radius of curvature is bounded away from zero on C, i.e., 

p 2 ^ )8 > 

on C. 

In [2, 5] the inequality 



(3.32) 



(3.33) 



?c(f)Wc' -1 



2ff d -f-cH 
ds ds 



M * +2 iAff* 



M^ffH^iKf)* (3 ' 34) 



was derived. Here t m denotes the minimum value of t on C, and c is given by (3.20). Inserting 
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(3.34) into (3.32) we then obtain (using 3.33) the inequality 



f BHxdy^'ZAilD-'b 



M(V)- 


d 2 V 


2 

ds 






'jr 



H 2 ^fH 2&+2 Hf) 2< 4 (3 - 35) 



Again, for a star-shaped boundary, we may takef i =x i {x l =x,x 2 =y) and e=0, since x 1 n i >0 is just 
the condition for star-shapedness. r 

It remains to determine a bound for H 2 dxdy in (3.12). Such a bound, however, was given 

JR 

in [1]; i.e., 



where k 2 is given by 



f H 2 dxdy ^k 2 j) H 2 ds = k 2 i V 2 ds 



k 2 =^ [4Xi (ff) m + c 2 + 2c inf { $ |«2Ai [ff] M * r]n . 

Ai 



(3.36) 



(3.37) 



The symbol ( ) M denotes the maximum value of the quantity in parentheses. 
Combining (3.35) and (3.36) and inserting into (3.12) we obtain the inequality 



I 



F 2 dxdy^2AA2D' 2 

R I J( 



M(V)-D(\-(t) 



d 2 V 
ds 2 






r 
c 

ds 



J ds ds ^^ 



ds 



ds 



+ - 



; H 



(3.38) 



The insertion of (3.38) into (3.5), (3.6) or the similar expressions for higher derivatives yields an 
expression which is not precisely in the form (3.4). However, (3.4) can be easily obtained if the 
arithmetic-geometric mean inequality is used in (3.38). This completes the proof of theorem I. 
In some cases it may be difficult to compute a particular solution <p and its derivatives in R and 
on C. In that case we may approximate u by any sufficiently smooth function i// and decompose 
as follows: 

u-i\f=V+W (3.39) 

where 

A 2 V=0 (3.40) 

in R and 

V=u-*lf 9 M(V)=M(u-ilf) (3.41) 

on C, while 

A 2 r = A 2 (w-i//) (3.42) 

in R and 

W=0,M{W)=0 (3.43) 

on C. 

We have already established a bound (3.4) for V^ p) (0). We seek now a bound for W {p) {0). 
Clearly, 

W(P)=\ G{P,P')tf{u-^)dxdy (344) 

Jr 



where G(P, P') is the Green's function in the simply supported plate problem, i.e., G(P, P') 
denotes the deflection of a simply supported plate at a point P' in R due to a unit point load at P. 4 
We now approximate W(P) by 



W(P)=\ [r(P,P')+£]A 2 U-i//)&dy, 
Jr 



(3.45) 



1 For a discussion of G cf. S. Timoshenko, Theory of plates and shells, pp. 156-167 (McGraw-Hill, 1940). 
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where 

r(P,P')=^ P ,logJp (3 .46) 

the singular part of G(/\ P'), and B denotes any biharmonic function (possibly zero). We approxi- 
mate apth derivative, W (p \P), by W P {P), where 



W P {P) 



[ Y{P,P')tf(u-xls)dxdy {P) + [ BP& 2 {u-x\f)dxdy. (3.47) 



Here [ ] (p) denotes the pth derivative with respect to the variables of the point P, and B p is an 
arbitrary biharmonic function, which may be taken identically zero. Clearly then at an interior 
point P 

W^{P)-W p {P) = [ [G{P,P')-T{P,P')Y p) ^{u-x\>)dxdy- [ B p A 2 (ii-i!))dxdy. (3.48) 

JR JR 

The order of differentiation and integration is interchangeable in the first integral since G(P, P') — 
T(P, P') is biharmonic in R. For simplicity we adopt the convention that for a function of the two 
points P and P', e.g., I"(P, P'), the symbol T {p) (P,P') is to denote differentiation with respect to 
the variables of the first argument. Thus T {p) (P, P') denotes differentiation with respect to the 
variables of the point P while Y {p) (P' , P) denotes differentiation with respect to the variables of the 
point P' . 

In view of (3.39) we obtain 

\i+ p \0)-^ p) {0)-W p (0)\ ^ |FW(0)| + \W<*KO)-W*(0)\. (3.49) 

The term W p (0) involves only known quantities and hence may be computed directly. If p ^ 3 
the differentiation indicated in the first term on the right of (3.47) may be taken under the integral 
sign and placed on r(P,P'). The resulting integral exists provided A 2 (tf — 1//) =0{r e P p, 1 ) as 
P—>P' for some positive e. If we wish to compute W p {0) for p ^3 we decompose W P (P) into 

F"(P)=F?(/>) + r$(P) (3.50) 

where 

W\ (P) = [ j R T(P, P')tf{u-*)dxdyY (3.51) 

and 

W\ (P) = J r fiPA 2 ( u - 4>)dxdy. (3.52) 

Since B p is biharmonic the last term presents no difficulty. In order to evaluate W?(P) it is neces- 
sary to make some assumptions on the smoothness of A 2 (u — \jj). We assume for the moment that 
A 2 (« — 1//) has a sufficient number of derivatives in some spherical neighborhood S s of radius 8 
about the origin and that P lies on the interior of Ss. Then 

W\(P) = jr^(P, P')A2(u-<l>)dxdy 

= | n i )(P,P')A"(B-*)dWy- f P»(P\P)A*(B-*)<fcdy. (3.53) 
Jr-s s Js s 

We now integrate by parts in S« and obtain 

W\{P)=\ r^(P,P')A*(u-<P)dxdy 
Jr-s & 

-l nT(P,P')&i(u-il,)ds+l r(P,P')[A2(u-4,)]^dxdy. (3.54) 
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Here we have used the fact that \(P, P') — Y{P' , P). In (3.54) Sg denotes the boundary of Sg and 
n' the component of the unit normal in the direction in which the differentiation is taken. Since 
P lies interior to Sg any derivatives of T(P, P') with respect to the variables of the point P are 
bounded on 2§ and in R—S?>. The last integral on the right of (3.54) may be differentiated 3 
times with respect to the variables of P and the order of integration and differentiation interchanged 
provided [A 2 (u — 1//)] (1) = 0(r)>p, l ) for some positive e as P' —> P in S§. Letting P^O we thus 
obtain the value_of W A {0) . 

To define W p (0) for p > 4 we differentiate (3.54), make use of the fact that r (,) (P, P') = — 
p(D(/>' ? p) ? integrate the resulting final term by parts, and continue in this manner until we obtain 
an integral over Ss which has the desired number of derivatives. We see in this way that for 
W p (0) to remain finite, it is sufficient that A 2 {u— i//) have p — 4 derivatives in Ss and that the deriva- 
tives of order p— 3 be 0{R e ~ p }) for some e > O as P' — > P in Sg. 

We need now only a bound for the last expression on the right of (3.49). To this end we use 
Schwarz's inequality in (3.48) to obtain 

[_W^(0)-Wp{0)Y^ [ {[G(0,P')-r(0,P')]W-BP}dxdy f [A 2 (u-xlj)] 2 dxdy. (3.55) 

JR JR 

Since the differentiations in (3.55) are with respect to the variable of P (before evaluation at O) 
the quantity {G(0, P') - T(0, P')]^-^} is just another function V, i.e., A 2 F= O in /^with Fand 
M(F) given on C. Thus (3.38) is valid with V replaced by { [G(0, P f )-T(0, P')] {1,) -Bp} through- 
out. Since G^(0,P') and M(G^(0,P f )) vanish on C only the quantity [T^(0, P') +5*>] 
appears for K on the right. We note that the right hand side of (3.55) may be made small by 
choosing B p to approximate the data of — T (p) (0,P') and/or by choosing i// to approximate the 
data of u. 

Let us denote the right hand side of (3.38) by I(V) and designate the constant in the mean value 
inequality for |FW(0)| as C p , i.e., 

1^(0)1 ^C p J Pd*dy. (3.56) 

Then our bound may be expressed in the simple form 

|^)(0)-^)(0)-^)(0)|^Cp/(a-^)+/(rp)(0,P , )+BP) f [A 2 U-*//)] 2 ^/y. (3.57) 

4. Third Boundary Value Problem for Elastic Plates 

In this section we determine a priori bounds for a function V, assuming A 2 V (the loading) is 
prescribed throughout a simply connected region R and that the boundary values of M(V) and (J(K) 
are prescribed on C. Here M(V) is equal to —M(V)/D where D is given by (2.1) and M(V) is given 
by (3.3). Also 



Q(V)=^(AV)-(l-cr) 
on 



dW d (\ bV 



ds 2 dn ds y> Bs 



(4.1) 



The quantities M(V) and Q(V), are proportional to the bending moment and the reaction. We 
prove in this section the following theorem. 

THEOREM II. Let V be any function with piecewise continuous fourth derivatives in R which 
satisfies the conditions 

c Jc ox Jc dy 
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Then at any point O in R, 

|FW(0)| 2 ^A;i<j> [M(V)] 2 ds + k 2 j> (Q(V)] 2 ds + fa f r)(A 2 V) 2 dxdy (4 2) 

where the constants k, are explicitly determined, p ^ 3, and r\ is a known nonnegative function in R. 

Here, as before, V^ p) (x,y) denotes any pth order derivative of V. As indicated in the previous 
section once an a priori inequality of the type (4.2) has been derived it may be used to obtain point- 
wise bounds in the indicated boundary value problem. 

Before proceeding further we define quantities p 2 and q 2 in terms of which the constants in 
(4.2) will be expressed. These quantities are defined as follows: 

p 2 = min - A — , (4.3) 

i^ =0 (j) x Hs 

where the minimum is taken over all piecewise continuously differentiable functions x which satisfy 
the condition (J) \ds = O; 

ds 

jc \os / 

qz= min — : (4.4) 



m 

i**- (f ^ds 



where the minimum is taken among all functions ifj continuously differentiable on C and satisfying 
<b \\)ds = 0. The eigenvalues corresponding to the Rayleigh quotient in (4.3) are commonly referred 

to as Stekloff eigenvalues. They are the eigenvalues of the following problem 

Ah = (4.5) 

in R and 

^-ph = (4.6) 

on 

on C. In general, p 2 will not be known explicitly. However, for our purposes any positive lower 
bound p 2 for p 2 will suffice. We indicate later how this lower bound is obtained. The eigenvalues 
corresponding to the Rayleigh quotient of (4.4) are the Wirtinger eigenvalues. They satisfy 

d 2 u , n 

^ qu = (4.7) 

on C and the proper continuity condition. The eq (4.7) is of course the equation of the vibrating 
string of length equal to the perimeter L of C. Thus 

_4tt 2 
Q2 ~W (4.8) 

We proceed now to the proof of theorem II. Let us again denote by T(P, P') the fundamental 
solution of the biharmonic equation. 



r(P,P')=^,iog(JL). (4 . 9) 
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(We could add any biharmonic function B to Y without altering the following results.) In terms of 
Y(P,P') we have 

V(P)=4> V^(T)ds-ij-Jl^)ds^ijl{V)j^ds-i ~Q{V)Yds+[ YtfVdxdy. (4 j ) 

For P on the interior or R we may differentiate (4.10) with respect to the variables of the point P 
and interchange order of differentiation and integration on the right. Thus we obtain for p ^ 3, 



V**\0)=i VQ{Y^)ds-i |^M(I»)ds + (p M{V) ^-Y^ds-b Q{V)Y^ds 
Jc Jc on Jc dn ] c 



+ 



f Y^&Vdxdy. 

JR 



(4.11) 



Here again as in the previous section we could add to T (p) any biharmonic function B p . The last 
three terms of (4.11) involve only the known function T and the data of V. Hence we define 

*fe>(0) = V&KO)-j> M(V)f^ds + j> Q(V)Y^ds- f Y^A 2 Vdxdy. (4.12) 
and obtain by Schwarz's inequality 



|(D(p)(0)| 2 



&v*ds+j> (j^fds {(^ r ) 2 f[(?(r (p) )] 2 ^ + (£ [M(r^)] 2 *j. (4.13) 



Since (p Vds = it follows from (4.4) and (4.8) that 



Thus 



Nov 



i V*ds ^ t 



mm-- 



|*<p>(0)| 



>(f) 



'Hm- 



£)% W(P p) )] 2 d s +j> [yw(P">)] 2 * 



,£>■*+ 



m-im + ffl. 



ds. 



(4.14) 
(4.15) 

(4.16) 



In view of (4.3) and the fact that CD I — J ds = d) ( — — J ds^O we have 



(4.17) 



Thus from (4.14), (4.16), (4.17) and the definition of E(V, V) we have for o- > 



m™ + m 



*Dp 2 (l-(j) 



E{V,V). 



(4.18) 



(For — V2 < a < O the factor (1 — a) is replaced by (1 + 2a).) From the divergence theorem we then 
obtain 

T)'£^ + £(© ,A *isn^{i^ sr(F) *-i^ (F) 4 (4 - i9) 
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An application of Schwarz's inequality yields 



m~+m** 



Dp2 a-<ri 



! {£ [M(V)Us+(^Jj> c [Q(V)ydsj- (4.20) 



Insertion of (4.20) back into (4.13) then gives the desired bound 
2 



|*<*>(0)| 



Dp 2 (l-<r) 



C [M(V)¥ds + &)^„[Q(V)¥ds} 



'{£ 



277/ JC 



{£ [M(I™)] 2 ds + ^Y£ [Q{rv>)¥ds\. (4.21) 



This is equivalent to the desired inequality (4.1), since instead of (4.12) we could have used 
Schwarz's inequality on the last three terms, being careful to use 

f Y^tfVdxdy^ \ r PP >[T^ydxdy f r? p A 2 Vdxdy for p = 3. 

JR JR JR 

It remains now to derive a bound for p 2 (i.e., find a p 2 ). We denote by u 2 the eigenfunction 
corresponding to p 2 , i.e., (see (4.4) and (4.5)) 



in R and 



on C. From (4.4) and (4.8) we have 



AlL 2 = 



du 2 n 

— -p 2 u 2 = 



i-^Yi l^\ 



(4.22) 
(4.23) 

(4.24) 



\2ttJ Jc\ds 
We consider again the identity (see [5]) 



which was used in deriving (3.34). Setting 

[f*k, 8 u~ 2/J i\ u 2 , [U 2 , idxdy ^ c u 2 , tu 2 , idxdy 

JR JR 

we obtain from (4.25) (by making use of the arithmetic-geometric mean inequality) 

* i&ff ds * i * ( **)* * * i (tf) 2 * + 2c L m ' m - idxdy ' 

where 

l = 2{ 



K /, aD ,+i, M * 2 ^ 26 *--'-> 



(4.26) 



(4.27) 



(4.28) 



with b given by (3.21). Inserting (4.28) into (4.24) and making use of (4.22) and (4.23) we obtain 

u 2 ds^ (4.29) 



^j\- m l {plb + 2cp,}^uHs. 
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It follows then that 



that 



p%+$ *4m*t m (bL*)- 1 +[cIb'}* 



(4.30) 



p 2 ^p 2 ^ [47r 2 t m (6L 2 )- 1 + (c/by]*-c/b. (4.31) 

If the boundary C of /? is star-shaped with respect to the origin (i.e., x l ni > 0) then we may 
choose/ 1 =x\ in which case c = 0, t = x i rii, and 6 = 2{2r 2 A/ £~ 1 — £ m }. 

Thus far we have indicated how to obtain bounds for the function and its first, second, and 
third derivatives. In physical problems this includes all derivatives for which we would normally 
desire bounds. However, if we wish bounds for higher derivatives they may be obtained using 
the techniques of the previous section. For interior bounds the only term which will cause diffi- 
culty for p > 3 is the last term on the right of (4.12); but this is precisely the term which was evalu- 
ated in the previous section. 
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